GLOBAL REGULARITY FOR ORDINARY DIFFERENTIAL 
OPERATORS WITH POLYNOMIAL COEFFICIENTS 

FABIO NICOLA AND LUIGI RODINO 

Abstract. For a class of ordinary differential operators P with polynomial co- 
efficients, we give a necessary and sufficient condition for P to be globally regular 
in R, i.e. u £ 6>'(R) and Pu g S(R) imply u G <5>(R) (this can be regarded as 
a global version of the Schwartz' hypoellipticity notion). The condition involves 
the asymptotic behaviour, at infinity, of the roots £ = £j(x) of the equation 
p(x,£) = 0, where p(x,£) is the (Weyl) symbol of P. 



1. Introduction and statement of the result 

Linear ordinary differential operators with polynomial coefficients play an impor- 
tant role in mathematics. On one hand, they were source for the study of classes 
of special functions, with links all around to Applied Sciences. On the other hand, 
the general study of Fuchs points, irregular singular points and Stokes phenomena 
present deep connections with different branches of geometry. 

We may write the generic operator in the form 

(1.1) P= a °,f> xf>D °'> a a ,peC, 

a+/3<m 

where we use the notation D = —id/dx. In this paper we deal with the Fourier 
analysis of (jl.ip . Namely, rather than analysing the extension of the solution 
u of the corresponding equation in the complex domain C, we shall address to 
their analysis in M? x ^ with respect to the phase-space variables x,£ G R. Basic 
function space in this order of ideas is the space «S(M) of L. Schwartz [17], defined 
by imposing 

(1.2) sup \x p D a u[x) | < oo, Vet, (3 G N. 

These functions are regular in M^, in the sense that both u(x) and its Fourier 
transform tt(£) = L e~ tx ^u(x) dx present a rapid decay at infinity, beside local 
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regularity. As universal set for our study we shall take the dual «S'(R); note that 
this will exclude solutions which have exponential growth at infinity. 

Aim of this paper is to characterize the operators P in (II. II) which are globally 
regular, according the following definition. 

Definition 1.1. We say that P is globally regular if 



In particular, if (II. 3ft is satisfied, the solutions u G <S'(R) of Pu = belong to 
«S(R). Global regularity turns out to be basic information in many applications. So 
for example in connection with Quantum Mechanics, assuming P in (II. ip is self- 
adjoint, we may deduce that the eigenf unctions, intended as solutions u G L 2 (R) 
of Pu = 0, are in <S(R). In the Theory of Signals, where we may regard P in (II .ip 
as a filter reproduced by electronic devices, we have that the globally regular P are 
exactly the ineffective filters, i.e. filters which do not cancel any essential part of 
the signal. 

The literature concerning global regularity, sometimes also called global hypoel- 
lipticity, in these last 30 years is extremely large, taking also into account the same 
problem for operators with smooth coefficients and pseudo-differential operators in 
R ra , with n > 1. We address to the recent monograph of the authors [13] for a 
survey. 

Let us begin with a simple example. The constant coefficient operator p(D) = 
J2a<m c &D a , c a G C, is globally regular if and only if p(£) ^ for all £ G R. In fact, 
if p(£o) = 0, then u(x) = exp[i^ x] provides a solution of p(D)u = 0, whereas by 
Fourier transform one gets easily that p(D)u G <S(R), u G <S'(R) imply u G <S(R) if 
p(0 7^ f° r an £ £ R- The same result keeps valid for partial differential operators 
with constant coefficients. 

Passing to operators with polynomial coefficients, a characterization of globally 
regular operators in R n is certainly out of reach at this moment. However for the 
ordinary differential operator (II. ip a necessary and sufficient condition seems pos- 
sible, and we shall give it in the following, under an additional algebraic condition. 

Consider first the standard left-symbol of P in (11.11) : 



(1.3) 



u G S'(R) and Pu G <S(R) 4n£ <S(R). 



(1.4) 




a+/3<m 
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In our approach, it will be convenient to argue on the Weyl symbol, see for example 
PH Chapter XVIII], given by 

(1.5) P(z,£)= c ^ x ^° 

a+/3<m 

We shall assume c m)0 = 1 in (11.51) . i.e. a m0 = 1 in (jl.4p . We have 

m 

(1.6) p(^o=n^-^)) 

where £_}(x), j = 1, . . . ,m, are real-analytic functions defined for x G M, |x| large 
enough. Let us denote by Ai, . . . , A m the complex roots of the polynomial 

m 

We have, possibly after relabeling, 

(1.7) £j(x)/x — > Xj as x G M, |x| — >• +oo. 

In fact, the roots have a Puiseaux expansion at infinity ([91 Lemma A. 1.3, page 
363]), namely 

(1.8) — XjX + Cj t k(x 1 ^ p ) k for |x| large, 

— oo<fc<p— 1 

for some integer p, where the function x 1//p is the positive pth. root of x for x > 
and, say, x 1 ^ = |a;| 1 ' /p e" r//p for a; < (by taking the the lowest common multiple 
we can assume that the same integer p occurs for every j). 

We suppose that the roots which approach the real axis at infinity are asymp- 
totically separated, in the sense that 

(1.9) If Xj = A fc G R, with j j£ k, then 

\€j( x ) - €k(x)\ > max{|^-(x) - Xjx\, \£ k (x) - X k x\, |x|~ 1+£ } 
for some e > 0. 

Theorem 1.2. Assume ( II .9p . Then P is globally regular, i.e. fll .3j) holds, if and 
only if 

(1.10) |xlm£j(x)| — > +oo as x G R, |x| — > +oo. 
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For a better understanding of fll.9p . ( ll.lOp . we may argue on the Puiseaux ex- 
pansions (11.81) . Let us emphasize the first non- vanishing term after Xjx, namely 
write 

(l.n) s./(-'') v • c^r^'y ■ Yl c ^( xl/p ) k 

— oo<fc<r(j) 

with Cj>(j) 7^ 0. Arguing for simplicity for x > 0, condition (11. 9ft states that, if 
Xj = Afc G M, then one at least between r(j) and r(k) is strictly larger than —1 and 
moreover, in the case r(j) = r(k) > —1, we have Cj jr u\ ^ c fcjr ( fc ). 

As for (ll.lOp . arguing again for x > 0, it means that for all j = 1, . . . ,m, we 
have ImAj ^ or else Imcj ^ 7^ for some k with k/p > —1. To be definite: (ll.lOp 
is not satisfied when for some j the function 

v+ c A xl/p ) k 
—p<k<p—i 

is real-valued. 

In the next Section 2, we recall some notation for the pseudodifferential calculus 
and we reduce the problem to the analysis of a global wave front set. The proof 
of Theorem 11.21 is given in Section 4, by using the factorization of the operator 
presented in Section 3. In some sense, this factorization is the counterpart in phase- 
space of the classical methods of the asymptotic integration. In fact, we think that, 
in terms of classical asymptotic analysis, under condition ( II. 9p . a proof of Theorem 
11.21 would be as well possible, but much more difficult. Section 5 is devoted to 
remarks and examples. In particular, we recapture some existing results on the 
global regularity for operators of the form (II. ip . giving corresponding references. 

We finally observe that Theorem 11.21 keeps valid for relevant classes of operators 
of the form ( II. ip . independently of the asymptotic separation ( II. 9p . For example, 
( II. 9p is not satisfied for constant coefficients operators with a double root = 
£k( x ) — £o> whereas in this class (ll.lUp is necessary and sufficient for the global 
regularity, as observed before. Regretfully, the assumption ( II .9p will be essential 
in our proof of Theorem 11.21 for operators with polynomial coefficients. 

2. MlCROLOCAL REDUCTION 

We recall that a pseudodifferential operator in 1R, according to the standard 
quantization, is an integral operator of the form 

p(x,D)u(x) = (2n)- 1 [ e^p(x,0md^ 

Jr 

for u G iS(lR), where the so-called symbol p{x,!;) is a smooth function in M? which 
satisfies suitable growth estimates at infinity, which will be detailed below. The 



5 



corresponding operator p(x, D) will define continuous maps <S(R) — > <S(R) and 
<S'(R) <S'(R). 

An important class of symbols is given by the space r m (R) of smooth functions 
p(x, £) in R 2 satisfying the estimates 



\d?dep(x,Z)\<C af} (i + \x\ + \Z\) m - a -f ) , Vc*,/3gN, 



r>2 



for some m G R. This class arises, in particular, in the following definition of the 
global wave- front set of a temperate distribution, as introduced in [TTj . 

A point (xo,£o) G R 2 \ {0} is called non- characteristic for p G T m (R), if there 
are e, C > such that 

(2.1) \p{x,Z)\>C{l + \x\ + \£\) m foT{x,Z)eV (xo ^ 

where, for z G R 2 , Zq ^ 0, V zo>£ is the conic neighborhood 

V zo . £ = {z G R 2 \ {0} : <e, \z\>e- 1 }. 

I |Z| |^o | J 

Let m G <S'(R). We define its (global) wave- front set WKF(u) C R 2 \ {0} by saying 
that (xo,^o) ^ R 2 ) 0eO)£o) 7^ (0,0), does not belong to WF(u) if there exists 
ip G r°(R) which is non-characteristic at (xq,£o), such that ip(x,D)u G «S(R). The 
set WF(u) is a conic closed subset of R 2 \ {0}. 

Proposition 2.1. ([H]) If the point (xo,£o) is non- characteristic for p, and (#o,£o) ^ 
WF(p(x,.D)u), ffcen (x ,fo) gWF(u). 

Proposition 2.2. ( [1 1J ) For every p G T m (R), andw G «S'(R), we haveW F(p(x, D)u) 
cWF(u). 

We will also use the following characterization of the Schwartz class. 
Proposition 2.3. ([TT]) For « G <S'(R) we /iai>e W / F(m) = if and only if u G 



In fact, results similar to those in Proposition 12.11 hold for more general classes 
of operators. The following case will be important in the following. 

Consider the class T™(R) of symbols p(x,£) satisfying the following estimates: 

(2.2) \d^ x p(x,0\ < C af3 (O m - a (x) m -P +Sa , V«,/3 G N, (*,£) G R 2 , 

for some m G R, < 5 < 1. For t G R, we write (t) = (1 + t 2 ) 1/2 . Notice that 
r°(R) C r°(R). These classes are a special case of the Weyl-Hormander classes 
S{M,g), with weight M(x,£) = (x) m (0 m and metric g x£ = fpr + {x) %& - The 
usual symbolic calculus with full asymptotic expansions works for these classes (cf. 
[T0| Chapter XVIII]), because the so-called Planck function h(x,£) = (C) -1 ^)' 5-1 
satisfies h(x, f ) < C(l + |ar| + l^) 5 " 1 and 5 - 1 < 0. 
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A symbol p G r™(M) is called hypoelliptic at (xq, £o) 7^ (0, 0) if for some e > it 
satisfies 

(2.3) \^ xP {x,0\<C a p\p(^0\(0~ a (^' P+Sa , Va,/3eN, (x,0eV M)>£ . 
and 

(2-4) \ P (x,C)\>C(xr'(O m ' V(*,O^Wo),e, 

for some m! G K., C > 0. 

Proposition 2.4. Letp{x,£) satisfy ffZTZty . fl2~31 and ([231). U{xq,^o) £ WF{p{x,D)u) 
then (x ,£ ) £ WF(u). 

The proof is standard, since the assumptions imply the existence of a microlocal 
parametrix Q = q(x,D), q G r^ m (IR) (see e.g. [51 Lemma 3.1] and [T31 Lemma 
1.1.4 and Proposition 1.1.6]). 

Finally, we will also make use of the the class r^ n (5 (IR) of smooth functions p G 
C°°(R 2 ) such that 

(2.5) \^p(x,0\ < C aP {l + \x\ + \i\) m - a+s ^ Va,/3 G N, (x,0 G K 2 , 
for some mG M, < 5 < 1. We have the following result. 

Proposition 2.5. Lei p G r^(R). T/ien ^F(p(i, D)u) C 

The proof is again standard and omitted for the sake of brevity. 
Finally we define, for future references, the following class of functions in R. We 
set S m (M.), m G M, for the the space of smooth functions satisfying the estimates 

(2.6) \f (a) (x)\ <C a {x} m - a } V«gN, xgR. 

As already observed in the Introduction, it is also useful to deal with the Weyl 
quantization, defined as 

p™(x,D) = (27T)- 1 J^-y^{^Y L ^)u{y)dydi. 

For example, an operator with a symbol a(x,£) as in (jl.4p in the standard quanti- 
zation can be always re- written as an operator with Weyl symbol p(x,£) given by 
the formula (11.51) . 

The main feature of Weyl quantization is its symplectic invariance: if x : M 2 — M 2 
is a linear symplectic map, there is a unitary operator U in L 2 (M) which is also 
an isomorphism of S(M.) into itself and of iS'(IR) into itself (in fact a metaplectic 
operator), such that 

(2.7) (po X ) w ( x ,D) = U- 1 p w (x,D)U. 



Moreover, in the definition of the global wave-front set it would be equivalent to 
use the Weyl quantization, which implies at once that 

(2.8) (x ,£ ) t WF(u) ^ x(x ,£o) € WF(Uu), 

where U is the operator associated to x as i n (12. 7p . 



3. Factorization of the operator 

In this section we provide a convenient microlocal factorization of the operator 
P in Theorem 11.21 which reduces the proof of that result to the case of first order 
operators. 

Let j G N and consider the elementary symmetric functions defined by 



o-q(£i, • 


••.&) 


= 1 




••,0) 


= - E 6 






i<i<i 


0-2(6, • ■ 




= E 






l<?<fc<7 


(7,(61,.. 


•■,&) 


= (-1)^1 ■■•« 



Proposition 3.1. Lei ri > 0, r 2 > and dj(x), j = 0, ...,r 1; and £j(x), j = 
r x , . . . ,r x + r 2 , fre smooth functions, defined in some open subset ofR. Then we 
have 

f>(*)LT- fc if (£>-&(*)) 

fc = J=F1 

n+ra 

= E( E ffl iW^(UW,-.^r 1+r2 W) + ^W)fl ri+r2 ^ 

fc=0 l+h=fc 

Z <ri M<.r2 

where Ro = R\ = and, /or k > 2, 

(3.1) i? fc G spaocia^j™ ■ • • ^ h) , Z + /i + mi + . . . + m h = k, 1 < h < k - 1, 

mi + . . . + m h > 0, r x + 1 < j x < . . . < j h < r x + r 2 }. 
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Proof. Let us apply induction on r 2 . The conclusion obviously holds with Rj = 
for every j if r 2 = 0. Suppose it holds for r 2 . The we have 



(3.2) E *j(x)D ri ~ k (D ~ Zr 1+ i(x)) ■ ■ ■ (D - ^ 1+ r 2+ i(x)) = 

fc=0 
ri+r 2 

E( E ^(xK(ei(a;),...,a(a;)+ J R fc (a;))^ +r2 - fc ( J D-e i+r2+1 (a;)), 



fc=0 i+h=fc 



with as in the statement. Let = ai(x)ah(£ri+i( x )i ■ ■ ■ i£n+r 2 { x ))- 

l+h=j 

By Leibniz' formula the expression in the right-hand side of (13 .2p reads 



n+r2+l 



(3.3) E ( & fc( X ) -Cri+r 2 +l(^)&fc-l(^) +i?fc(x)- 

E ^ + r2 + 1 ~~ * " J ) (bk-j-iix) + R k ^ l (x))D^ ri+r2+1 (x))D^+ 1 



3=1 

Now, we have 



h(x) - £ ri+r2+ i(x)b k -i(x) = E ai(x)a h (£ ri+1 (x) } . . . , £ ri+r2+1 (x)) , 



l+h=k 
I<r 1 ,h<r 2 + 1 



whereas the terms Rk(x), bk-j-i(x)D 3 £ n+r2+1 (x) and i?fc_ J _i(a;)Z) J ^ r . 1+r2+1 (x) are 
admissible errors. □ 

Proposition 3.2. Let r\ > ; r 2 > and Oj(:r) ; j = 0, ...,r\, and bh{x), h = 
0, . . . , r 2 6e smooth functions on M. Consider the operator P with Weyl symbol 

(3.4) P (x,o = E%(x)r^' 

i=o /i=o 
Then P in the standard quantization has symbol 

qM= E { Rk+ E a ^) e 



fc=0 ' i+h=k 



Rk G span c {(aj&/i)^, 1 < i> < k, I + h = k — v\. 
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Proof. From (pQ| we have p(x^) = YJk=Q 2 J2i+h=k a i b h C 1+r2 ' k - On the other 
hand, for the standard symbol g(x, £) we have the formula 



u=0 

so that the desired result follows at once. □ 
Lemma 3.3. Let ^GC, k — 1, . . . , r\ + r<i, 

Qio = n ^ - & 

ri+l<fc<ri+r2 

and 

q,(o - n 

l<fc<r^+r2 

for j = n + 1, . . . ,ri + r 2 . 

Consider the square matrix A of size r\ + r 2 , whose jth column for 1 < j < r\ is 
made of the coefficients of the polynomials £ ri_ -'Q(£)? whereas the jth column for 
7*1 + 1 < j < ri + r2 made o/ iae coefficients of the polynomial Qj (hence in the 
jth row there are the coefficients of x ri+r ' 2 ~ : ' which appear in those polynomials). 
Assume that ^ £j when j ^ k. Then the matrix A is invertible and its inverse 
A^ 1 = B = (Bjk) is given by the following formula: 

r l c n+r 2 -k 



(3.5) B jk = \ — — — — — (7j-i(^i, . . . , £a_i, ^+i, ■ ■ ■ , £ n 

l<Z<r^+r2 



if 1 < j < r l5 1 < k < r\ + r 2 , whereas 

*r\+T2—k 

(3.6) £ 



l<i<r 1 +r 2 

if ri + 1 < j < ri + r2, 1 < k < r\ + r%. 

Proof. To compute the fcth column of the inverse matrix, we have to solve the 
system AX = Y, where Y = [0, . . . , 0, 1, . . . , 0]* is the feth element of the canonical 
basis of C ri+r2_1 . Now, the vector X = [ai, . . . , a n+r2 ]* is a solution if and only if 

(3.7) r(oq(o+ E s^(0 = r i+r2 - fc , 

j=n+l 

with i2(e) = E;Li%e n - j '- 
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To solve (13. 7j) . we express the right-hand side in terms of the interpolating La- 
grange polynomials associated to the points j = 1, . . . , t*i + r 2 , which are 

L h (0-= II /i = l,...,r 1 + r 2 

1 <fc<r-j^ +r2 

(hence L/i(a*) = Qh(x) for /i = r 1 + 1, . . . , r 1 + r 2 ). We have 

(3.8) r+ r 2 - fc= £ ( ^ 

hence we are reduced to solve (13. 7\\ when the right-hand side is replaced by each 
Lagrange polynomial. On the other hand, one sees immediately that the equation 

ri _i_ r2 

R(0Q(0+ E a J Q J (0 = L h (x), 

j=r i+l 

when 1 < h < r 1 has the solution R(£) = Yl (6 — 6)> an d «j = for every 

i<;<r x 

j — n + 1, . . . , ri + 7*2, whereas when ri + 1 < /i < ri + r 2 the solution is given by 
R(x) = and = 1, Oj = if 7*i + 1 < j < r\ + r 2 with j ^ h. 

By linearity and (I3.8P we obtain that the solution of (13 ,7p is given by 

[8(8 = ^,^1 K-6) 

£ r l+ r 2- fc 

I/O "/., • J =r 1 + l,...,r 1 +r 2 . 

We can now come back to the original system AX = Y, and we obtain the desired 
form for the inverse matrix. □ 

It follows from the proof of the following proposition that the result in Lemma 
13.31 continues to hold if one just require that £ 3 - ^ for the indices j ^ k such that 
j > n or k > r\. In fact, as we will see, a simplification of the factors ^ — with 
1 < k 7^ / < r\ occurs in (13. 6p . 

Lemma 3.4. Let £j(x), j = 1, ■ ■ ■ , r± + r 2 fre smooth functions for x large, such 
that £j(x) ^ £k{%) tfn + 1 < j ^ k < n + r 2 , 

(3.9) Zj(x) = 0{x), if 1 < j < n + r 2 
and 

(3.10) ICj'C^) — CftO 3 ?)! ^ x if 1 < i < 7*i and 7*i + 1 < < 7*1 + r 2 
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as x — > +00. If we set £j = £j(x) in Lemma \3.3i for the inverse matrix Bj k 
Bj k (x) given there, the following asymptotic formulae hold as x — > +00: 

(3.11) B jk (x) = 0(xi- k ) 

for 1 < j < r\, 1 < k < r\ + r-i, whereas 

u*)\ ri+r2 - k 



(3.12) B jk {x) = 0\ 



n Mx) - &X) 

r-^+l<Z<r-^+r2 



for r% + 1 < j < r% + r 2; 1 < k < r% + r 2 . 

Proof. Consider first the entries with 1 < j '• < r%. With the notation of Lemma [3~31 
we can rewrite (13. 5p as 
(3.13) r 

Ef-iinr^-ife,..,^,^,.,^) n Q&) n (&-&: 

h=l i<;< ri i<i<i'< ri 

„ l,l'yth 

jk = n &-&o n om 

l<Z<«'<n l<l<n 

We claim that the numerator of this fraction, as a polynomial in £1, . . . , £n+r 2 is 

divisible by the product f] (£j — £;/) which arises in the denominator. To 

l<W<n 

this end, it is sufficient to show that the numerator vanishes if £ M = £„, for every 
1 < // < ^ < r\. It is clear that all the terms in the above sum with h 7^ /1, v 
vanish if £ M = £ u , because of the presence of the factor £ M — £ u . Let us prove that 
also the sum of the two terms corresponding to h = fi and h = v vanishes. Due to 
the symmetry of <x,_i, it suffices to show that 

(-ir 1 II (6-60 + (-ir _1 II (6-60 = o. 

l<;<;'<rj l<i<i'<r 1 
i.i'^M 1,1' 

This amounts to prove that 

i<i<n i<i<n 

where we took into account that fi < v. But this is true, because the two products 
that arise in the last formula coincide when ^ = £„. 

This proves the claim and shows, by a limiting argument, that the the matrix A 
in Lemma [3731 is still invertible if one just require £j 7^ for the indices j 7^ k such 
that j > ri or A; > ri. 
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Now, by dividing the numerator of (13 . 1 3f) by Yl (£t — 60 we S e t a polynomial 

i</<i'<n 

in £l, • • • ,fri+r 2 of degree^ 

(n + r 2 - fc) + (j - 1) + (n - l)r 2 + (ri ~ - 2) - ri(ri 2 -1) = r x r 2 -k+j. 

Hence, when we compose such a polynomial with the functions £ 3 - = by ( 13 .9p 

we get a function which is 0(x Tir2 ~ k+: >) as x — » +oo. On the other hand, using 
( EE} we have | J] Q(^i)| £ % rir2 , so that we deduce (ETTT1) . 

l<Kn 

Finally, let us prove (I3.12p . This follows at once from H3 .61) taking into account 
that, if ri + 1 < j < r\ + r 2 , 

n i^)- >^ n 

l<i<r^+?'2 r^+l<Z<rj +r2 

which is a consequence of ( 13. 10ft . □ 

Proposition 3.5. Letp(x,£) be as in ( 11.51) . Let £i(x), . . . , £ ri+r2 (x) ; ri + r 2 = m, 
be its roots, defined for x large, and assume (I3.9p . (I3.10p . as well as 

(3.14) >max{|^-(a;)|,|^)U" 1+e }, for j, = n + 1, . . . , n + r 2 , 

/or some £ > ; as x — > +oo. Let P be the operator with Weyl symbol p(x,£). We 
can write, for x large, 

ri n+r2 

(3.15) P = J2vk(x)D r '- k (D-^ar^ + fl, 

fc=o i=-ri+i 

where r\ k e S J '(K), = 1, . . . ,r 1; 60^) = 1; % £ ^(M), j = ri + 1, . . . T\ + r 2 (see 

(3.16) Vj (x) = *jfa(x), . . . ,£ n (x)) + 0(x>- 2 ) for j = 1, . . . ,n 
and 

(3.17) = + ©(a; -1 ) for j = n + 1, . . . , T\ + r 2 

as x +oo ; whereas R is a differential operator of order r± +r 2 whose coefficients 
are in S m (M) and rapidly decreasing as x — > +oo ; together with their derivatives. 



1 Observe that the polynomials J J (£; — and J J (£/ — £;/) have degree 



(r 1 -l)(r 1 -2) 



i<i<i'< ri i<l<l'<n 

l,l' 



and rii,r \ ^ respectively. 
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Proof. The standard symbol of P is given in Proposition 13.21 On the other hand, 
the symbol of the first operator in the right-hand side of (13. 15j) is computed in 
Proposition 13 . 1 1 By equating the coefficients of the terms of the same order we 
deduce that, modulo rapidly decreasing functions as x — > +oo, it must be 
(3.18) 

r\+r 2 

Y Y a l a h(£rx+l> ■ ■ ■ > £r 1+ r 2 ) + R'k = Y Vl( X )°h(Vri+l, ■ ■ ■ Vr 1+ r 2 ) + K 

k=0 l+h=k l+h=k 

for k = 1, ... , ri+r 2 , where the 77/s are unknown, a t = ai(x) = ai(£i(x), . . . , £ ri (x)), 
R' k e span c { (aia h (C ri+ i, £ ri+r2 )) M , 1 < v < k, I + h = k - u}, 

and 

Rl E span c {a^ ( ™ l} • • • r/j™' l) , I + h + m x + . . . + m h = k, 1 < h < k - 1, 

m 1 + . . . + m h > 0, ri + 1 < j x < . . . < j h < r l + r 2 }. 

We set 



(3.19) 



r]j{x) = aj(x) + (j(x) for j = 1, . . . , r x 

Vj( x ) = £j{ x ) ~ Gj(x) for j = n + 1, . . . , ri + r 2 , 



in the right-hand side of (I3.18p . By isolating the terms which are of order or 1 
with respect to ( = . . . , ( ri +r 2 ) we can write 

Y Vl Cr h(Vr 1 +l,---Vr 1 +r 2 )= Y a l^h(^r 1+ l, ■ ■ ■ ,^r 1+ r 2 ) + (AQ k + X k , 
l+h=k l+h=k 

where A is exactly the matrix described in Lemma [3. 3 1 By multi-linearity we can 
also write R'l = Yk + Zk, where Zk is the term constant with respect to the Q's. 
The system (13. 18[) hence becomes 

(3.20) C = -A^X - A~ l Y - A~ X Z - A^R'. 

Let us prove that 



(3.21) (A-'Z) 
as x —7- +00. We have 



j 



0(xJ- 2 ) ifj = l,...,n 

0(x~ l ) if j = ri + 1, . . . , 7*1 + r 2 . 



(3.22) Z k E spanja^j™ ■ • ■ $™ h \ I + h + m x + . . . + m h = k, 1 < h < k - 1, 

m 1 + . . . + m h > 0, ri + 1 < j x < . . . < j h < r x + r 2 }. 
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Hence, if j < n and A" 1 = (B jk ), we have (A~ 1 Z) j = EILT' B jkZ k , and by (13T221) . 
flU} and (IXTTD 

because k — I — h > 1. 

If j > n + 1, since = 0{^{x)x- m ») and m x + . . . + m h = k - k - 1, by (l3~22l) 
it suffices to prove that 

1 ' ) x ri n 

r-j^ + l<i <r-^ 

Mi 

By (J3H we have ^(x) = 0(\£ jv - if 3v + J and = 0(\^(x) - ^{x)\) 

and x\£j(x) — £,i(x)\ — >■ +oo if j ^ Z, i > ri + 1- Hence we obtain 

(3.24) x -n +k - h -i . Ux) ^~ k ^ X }'''^[f = 0(1) 

Mi 

(this is easily verified by considering separately the cases when the number of 
factors in the numerator is less or greater than that of the denominator). On the 
other hand, by Q we have \ ai (x)\x- k+k+l x- k+h+1 = 0(x~ 2ik - l - h)+1 ) = O^" 1 ), 
which combined with (I3.24p gives (13.231) . This proves (13.211) . 
Similarly one proves that 



(3.25) (A^R" 



0( X J- 2 ) ifj = l,...,n 

0(x~ 1 ) if j — ri + 1, . . . , ri + r 2 



as x +oo. 

Now, (13.211) and (I3.25f) suggest to look for functions Q with asymptotic expansion 

+oo 

(3.26) Q(x) ~ ^c J> ar'~ 2 ~p, for j = l,...,r u 

(3.27) Q(x) ~ J^c^a;" 1 "^, for j = n + 1, . . . , n + r 2 , 

where p is the integer for which (II. 8p hold. In order for the system (13.201) to be 
solvable by an iterative argument (in the sense of formal power series), it suffices to 

prove that the coefficient in (A~ 1 X)j, respectively (A~ 1 Y)j, of ar' -2- ?, respectively 

x p , depends only on the c JiAt with jj < k. To this end, consider first the term 
(A~ l X)j. By the very definition of X k , it belongs to the complex span of 
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with l + h = k, I < n, h < r 2 , where 5 t = a t or 5 t = Q, and 7^ = or 7^ = 
and the above product contains at least two factors of type (. 

Let 1 < j < t x . We have {A~ l X) J = Ek=P B jkX k . By (J3HJ, ([321 
and f)3.27p we get Bj k X k = 0(x^~ k x l+h ~ i ) = 0(x^~ A ). More generally, the same 
argument shows that for a fixed /i > the coefficients Cj )At , i = 1, . . . , ri + r 2 , 
may appear in the asymptotic expansion of (A~ l X)j only in the terms of degree 
< j _ 2 - H - 2. 

Let now 7*1 + 1 < j < 77 + r 2 . Let us prove Bj k X k = o(a; _1 ). It suffices to prove 
that 

(3.28) WfoWI^-SW-Ti.W = 

i ri n i&m-smi 



if I + h = k, I < ri, h < r 2 . We can suppose that, say, 7^ = Q h . Assume 
furthermore that there exists 1 < v < h — 1 such that 7^ = Then we have 
5i{x)x^ k+h ^jj h (x) = 0(x~ l ) and 



n u*)- 6(^1 0(1) ' 



7'-^ + l<Z<r-^ + 



where we used the same arguments as in the proof of (13.241) and the fact that 
Q v (x)/Mx)-^(x)\ = O(x-Vl0(^)-^(^)i) = o(l). Thisgives (ESSD- If instead 
there is no such a z/, it must be Si(x) = (i(x). Hence 5i(x)x~ k+h+1 = 0(x~ l ) and 

Ux)\ r ^~ k ln (x) ■ ■■ ljh (x)x-^-h-i 



n ux)-^x)\ 



which still gives (I3.28p . This proves that, if ri + 1 < j < r\ + r 2 , (A~ l X)j = o(x~ l ) 
as x — > +00. More generally, the same argument shows that for a fixed /z > 
the coefficients i = 1, . . . , ri + r 2 , may appear in the asymptotic expansion of 
(A~ 1 X)j only in the terms of degree < — 1 — -. Similary one sees that the same is 
true for (A~ 1 Y) j . 

Hence, the system (I3.20p has a formal power series solutions in the form ( I3.26p . 
f )3.27p . as x +00. Now, by a classical Borel-type argument, see e.g. [TUl Propo- 
sition 18.1.3], one can construct functions Q G S- 7 ' -2 ^) if 1 < j ' < 77, (j G S' -1 (R) 
if ri + 1 < j < r\ + r 2 , with the asymptotic expansions in f!3.26j) and (I3.27P 
respectively. 

The functions rjj (x) defined in (I3.19P for x large and extended smoothly to zero for 
x < then fulfill ( 13 . 1 5[) for a convenient operator R having the desired properties. 
This concludes the proof. □ 
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4. Proof of the main result (Theorem 11.21 ) 

4.1. Sufficient Condition. Let u be in S'(M) with Pu G <S(R). Let us prove 
that, under the assumptions (jl.9p and fll.lOp we have u G <S(R), i.e. WF(u) = 
(Proposition I2.3j) . 

We use the factorization (11.61) for the Weyl symbol p(x, £), valid for large \x\, and 
define the complex constant Xj, j = 1, ... ,m by (II. 7p . Let now (x ,^o) (0,0). 
If (x ,£o) does not belong to any of the rays {t(l, Xj), t G R \ {0}}, with Aj G R, 
we see from (11.61) . (I1.7P that p is non-characteristic at (x ,^q), i.e. it satisfies the 
estimate (12. ip for some e > 0. Hence (xq,^q) ^ V^F(u) by Proposition 12.11 

Let now (xo,£o) ^ e on a ra Y ^i)> £ G R\ {0}}, for some Aj G R. We can 
suppose that (xo,£o) = 0-7 ^j) or ( x o?£o) = ( — 1, — Xj). We can further reduce to 
the case when (xo, £o) = (1, 0). In fact, suppose that (xo, £o) = (1, Xj), and consider 
the linear symplectic transformation x(x,£) = (x, £ + Ajx). Let U be the operator 
associated to \ via (12 -7p . By (12. 8p . in order to get (xo,£o) ^ VKF(ii) it suffices 
to prove that (1,0) £ WF^U^u). Now, we have (p o x) w (x, D)U~ l u = U^Pu G 
«S(R), so that our original problem is equivalent to a similar one with (x ,£o) 
replaced by (1, 0) and the symbol p(x, £) replaced by (p o x)(x, £). The same holds 
for (xo,£o) = (— 1, — A) if we perform the preliminar symplectic transformation 
x(x,0 = (-x,-£). 

After these transformations we get a symbol, which we will continue to call 
p(x, £), which has a factorization as in (I1.6P for x large, where the new roots £j(x) 
satisfy (II. 7p and (I1.8P for other values of Xj, as well as (I1.9P and f ll.lOp . 

Suppose now Xj ^ for j = 1, ... , r\, and Xj — for j — r% + 1, . . . , T\ + r 2 , with 
r i + r 2 = Accordingly, by (II .7p . (II .9p we have 

(4.1) | x x for 1 < j < ri, 10(^)1 = o^) f° r Ti + 1 < J < Ti + r 2 , 

(4.2) |e,-(x) - £ k (x)\ > max{|e i (z)|, ^(x)!, x~ 1+e } for n + 1 < j < n + r 2 , 

as x ->■ +oo. Let us verify that (1, 0) £ WF(u) if Pw G «S(R). 

We can apply Proposition 13.51 and use the factorization in (I3.15p . We claim 
that the symbol of the first factor, namely Y7k=o^ k ( x )^ ri ~ k ; > * s n on-characteristic 
at (1, 0). In fact by (14. ip we have 

jn ri 

\j2^i( x )7---,^m n ' k \ =111^-^)1 for e v ^ 

k=0 j=l 

if e is small enough. Now by (I3.16P we deduce that 

ri 

IX)^)^ 1 "*!^ 1 " 1 ' for (x, o e y (1 , 0)l£ 

k=0 
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which proves the claim. Hence, since Pu G S(M) and (1, 0) G" WF(Ru), by Propo- 
sition |5TT] and A3. 15 j) we deduce that (1, 0) does not belong to the wave-front set of 
rifc r [+i(^ — r]j(x))u. Hence, to finish the proof it is sufficient to prove that, for 
every j = r x + 1, . . . , r x + r 2 , if (1, 0) g WF((D - ■q j {x))u) then (1, 0) WF{u). 
Fix such a j. We have from (jl.8p that 

(4.3) £j(x) = c j)k x k/p for x large. 

— oo<fc<p— 1 

By (13.171) we have an asymptotic expansion 

(4.4) r]j(x) = Cj^ k x k/p + 0(x _1 ) for x large. 

— p<k<p— 1 

By (ll.lOp . there exists > — p such that c^js. ^ R. Let < p — 1 be the greatest 
index for which this holds. Let Q(x) be a smooth real- valued function, Q(x) = 
for x < 0, and Qj(x) = f Y2v<k<p-i 

c jtk t k/p dt for x large. We can write 

(4.5) D - nj(x) = e iQ ^ x) {D - %(x)) o e ' iQ ^ x \ 
where % G 5 1 (R), and 

(4.6) %(x)=^(x)- c J,kX k/p = c hV x vlv + o[x vlv ) 

v<k<p—l 

for large x. 

Now, we can regard e^^i\ x ) as symbols (independent of £) in the class r ( j ) <5 (R) 
with 5 = 1- l/p (see ( 12. 5ft ). Hence, if (1,0) G" WF((.D - r]j(x))u), by Proposition 
12.51 we get (1,0) G" W^F^-D — fjj(x))u). On the other hand, the symbol £ — %(x) 
belongs to fJ(R) (see (E2D), and it is hypoelliptic at (1, 0) (see (IZSjl .lCT). In fact, 
using (14 .6p and the fact that Cj )U G" R we have 

\d?d^-Vj(x))\ <C Q/3 |£-%(x)|(£)- a (x)-^, Va,/3GN, x > 1, £ G R, 

with S = max{0, — v jp\ < 1 — l/p < 1, and |£ — %(x)| > |£| + x"/ p for x large and 
£ G R. Hence, by Proposition 12.41 we obtain (1,0) G" WF(u), which concludes the 
proof. 

4.2. Necessary Condition. Let us assume (II .9p but suppose that (ll.lOp fails 
for some j. Let us prove that then there exists u G <S'(R), u G" <S(R) such that 
Pu G <S(R). By the applying the same arguments as in the proof of the sufficient 
condition we can assume that (ll.lOp fails, say, for x — > +oo and j = r\ + r 2 , with 
X ri + r2 = 0. Hence (14. 3 p and therefore (14. 4p hold with j = r\ + r 2 and c ri+r2 ^ G R 
for every — p < k < p — 1. Let us set 

Q{x) = / ri ri +r 2 (t) dt, 
Jo 
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and let u G C°°(R) such that u(x) = for x < 0, u(x) = e %c ^ x > for x > 1. Then, 
by f)4.4p with j — r± + r 2 we have x~ c < |^(^)| ^ £ c when x > 1, for some c > 1. 
Hence w G <S'(R), m G' <S(R). Moreover, Pu = for a; < 0, whereas for x large 
we have (D — r] r . 1+r2 (x))u = and therefore, from (13.151) . Pu = Ru. Now, we also 
have \d a u(x)\ < x c+ ^ when x > 1, because ^n+r 2 G 5' 1 (R), so that is rapidly 
decreasing as x — >• +oo together with its derivatives. This implies that Pu G <S(R) 
and concludes the proof. 

5. Examples and remarks 

We give in the following some examples and remarks. 

Remark 5.1. The assumption c m o for the Weyl symbol of P can be eliminated, 
provided we submit preliminarily P to the conjugation by a metaplectic operator 
U. Namely, let U be associated to the symplectic map x( x i£) = ( x ~ £)■> A G R. 
In view of (12 .7p . the Weyl symbol of P = U~ l PU is given by 

(p°x)0&,o= E ^(x + ao^ 

o+/3<m 

with new coefficient 

Cm,0 = 2J 7^ 

Q+/3=m 

for a generic choice of A. On the other hand, the global regularity of P and P 
are equivalent, because U, U' 1 : S(R) -»■ <S(R), C/, f/" 1 : «S'(R) <S'(R). We may 
correspondingly re-formulate condition (I1.9P and Theorem 11.21 

Example 5.2. (Globally elliptic operators, cf. Grushin [5], Shubin [T5], Helffer 
[7J). Assume all the roots Xj, j = l,...,m, in (11.51) satisfy ImAj ^ 0. This is 
equivalent to the so-called global ellipticity of the symbol: 

(5.1) |p(z,OI >e(l + |z| + |£|) m for |x| + \t\ > R, 

for some e > 0, R > 0. In fact, if we write p(x, £) as a sum of homogeneous terms 

(5.2) p(x,£)= Pi( x >®> PiM= E c ^^, 
the condition (15. ip is equivalent to 

(5.3) p m (x,Z)^0 for (x, 0^(0,0). 

The corresponding operators are globally regular, since (ll.lOp is obviously satisfied. 
Basic example is the harmonic oscillator of Quantum Mechanics: 

P = D 2 + x 2 . 



19 



Example 5.3. Consider now the case when Xj G R for some j. Assume for the 
moment that all the real roots Xj are distinct, that is d^p m (x, Xjx) 7^ for x 7^ 0. 
We may apply Theorem II .21 By factorization we have 

£j(x) = Xjx + Cj + 0(x~ 1 ) 

with Cj = p m _i(x, Xjx)/d^p m (x, Xjx). Then P is globally regular if and only if 
Imcj 7^ for every real Xj. Consider for example the elementary polynomial 

P(x,0 = f - x + c. 

Limiting attention to the corresponding homogeneous equation, a classical solution 
of Du — xu + cu = is given by u(x) = exp[ix 2 /2 — icx}. If Imc 7^ 0, then u (jL «5'(M), 
whereas if c G R, then u G «S'(R), u <S(R), contradicting the global regularity. 

Example 5.4. (Quasi-elliptic operators, cf. Grushin j6], Boggiatto, Buzano, 
Rodino [2j, Cappiello, Gramchev, Rodino [4]). We pass now to consider the case 
when two, or more, real roots coincide. For simplicity we shall assume Xj = for 
all j — 1, ... , m, that is c a ^ = for all a, (3 with a + (3 = m, apart from c mi0 — 1- 
We may then consider the largest rational number q > 1 such that a + q/3 < m for 
all (a,/9) with c a ^ 7^ 0, and write 

(5-4) P(x,0= E c ^r, 

a+g/3<m 

with Cq,^ for some (ct,/3) with a + g/3 = m, a 7^ m. For the moment, we 
understand q < oo, that is a term with / actually exists. Note also that 
1 > ~~ n h\- Ths symbol in (15. 4ft is called (globally) quasi-elliptic if 

(5.5) Pm,,(a:,0= ^ c «/^^° for (x, ^ (0, 0). 

This implies in particular c 0>m / g 7^ and m/q is a positive integer. The correspond- 
ing operators are globally regular, as proved for example in [6], [2], |4]. Computing 
the Puiseaux expansion (11.81) . we have 

(5.6) ij{x)=rf\x\ l/q + o(\x\ 1/q ) for x^ ±00. 

where are the roots in C of p mg (±l,r) = 0. We deduce that (15. 5p is satisfied if 
and only if Im r- 7^ for all the roots. 

We then recapture the global regularity from Theorem ll.2[ provided condition 
(II. 9p is satisfied, i.e. the roots r~j~, or equivalently rj, are distinct, cf. (II. lip . As 
example, consider the Airy-type operator 

p = D 2 + cx, 

which is globally regular if and only if Imc 7^ 0. 
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Example 5.5. Let the symbol p{x, £) be of the form (15.41) . with distinct roots 
in C I5.6j) . but assume now r+ G R, or r~ G R, for some j. The condition (II. 9ft 
is satisfied, and we are led to determine the subsequent terms in the Puiseaux 
expansion (11. lip . We refer to Bliss [1] for general rules of computations, and we 
limit ourselves here to the example 

(5.7) P (x,o = r-^r-x 

with < r < m, A G R. The expansion in (I5.6P reads in this case 

£i(x) = e jX 1/m + o(\x\ 1/m ) 

where ej, j = 1, . . . , m, are the mth roots of 1. To fix ideas, assume m even; then 
we have for x — > +oo the two real roots ±1. It is easy to compute 

f ± (£) = ±x 1/m + c±x s + o(x s ), x — y +oo, 

with s = (r + 1 — m)/m and c± = ±A/m. Note that 1/m > s > —1. Hence if 
ImA ^ 0, the condition (ll.lOp is satisfied and Theorem 11.21 gives global regularity. 
Let us test this result on the corresponding homogeneous equation 

(5.8) D m u - AD r u — xu — 0. 

Every solution u G <S'(R), or u G <S(R) of the equation can be regarded as inverse 
Fourier transform of a solution v G «S'(R), respectively v G «S(R), of 

Dv + (x m - Ax r )v = 0, 

having the (classical) solution 

v ( x j = e i J 4 a; '-+ 1 /(r+l)-« m + 1 /(m+l) 

If 1mA ytz o, then v G" <S'(R) or v G «S(R), depending on r and A, that agrees 
with the global regularity of the operator. If A G R, then v G «S'(R), v £ <S(R), 
contradicting global regularity. 

Example 5.6. (Multi-quasi-elliptic operators, cf. Boggiatto, Buzano, Rodino 
[2]). By conjugation with Fourier transform, which we may consider as a meta- 
plectic operator, we can treat operators with symbol of the form (15 .4[) where the 
role of x and £ is exchanged. Global regularity is granted by ( 15. 5 \ or (15. 6p with 
Im ^ where we exchange again x with £ ; relevant examples of the correspond- 
ing operators are 

D + ix m , D 2 + x 2m , 

for m > 1. Multi-quasi-elliptic symbols are products of the symbols of this form, the 
quasi-elliptic symbols in Example 15.41 and the globally elliptic symbols in Example 
15. 2| possibly perturbed by terms in the interior of the Newton polygon generated 
in this way, see [2] for details and equivalent definitions. Under the condition ( 11. 9p . 
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we recapture from Theorem II .21 the result of global regularity in [2j. Limiting again 
to an example, consider the symbol 

£ 3 + ix£ 2 + x 2 . 

We have £1(2;) = -ix+o(\x\), £ 2 (x) = Vix 1/2 +o(\x\ 1/2 ), £ 3 (x) = -Vix 1/2 +o(\x\ 1/2 ), 
writing ±\/i for the two roots of i. Hence (ll.lOp is satisfied, and the corresponding 
operator is globally regular. 

Example 5.7. (SG-elliptic operators, cf. Parenti [T4], Cordes [5], Schrohe [TB] . 
Schulze [TH Section 1.4]). The case q = 00 in Example 15.41 corresponds to the 
case of the operators with constant coefficients, that we treated in Section 1. More 
generally, we can consider symbols of the following form, with m > 0, n > 0: 

(5.9) p(x,£)= ^Tc^C 

a <m, 
/3<n 

where we assume c m ^ n = 1. We say that the symbol ( 15. 9p is SG— elliptic if 

(5.10) \p(x,Z)\>e(l + \x\) n {l + \Z\) m , for \x\ + \£\ > R, 

for some e > 0, R > 0. In fact ( I5.10p is equivalent to the following couple of 
conditions: 

(5.11) ^aVO for^O, 

(5.12) J2 C ^ X ^° fora; ^ - 

/3<n 

We know from [HI [5j [151 HH] that SG-elliptic operators are globally regular. Willing 
to apply our Theorem II .2\ we first observe that factorizing p(x, £) we obtain the 
roots 

(5.13) €j(x)=rj + o(i), j = l,...,m, 

where r 3 - are exactly the roots of (15. lip . Since ( 15.1 ip is equivalent to Im rj ^ 
for every j, (15. lip implies fll.lOp for ( I5.13p . Similarly we can argue on the local 
ellipticity condition ( I5.12p . by using Remark IBTTl So we recapture global regularity, 
in the case when (II. 9p is satisfied (i.e. the roots of ( 15. lip . (I5.12p are distinct). 

Example 5.8. When, fovp(x, £) as in ( I5.9p . the equations (15. lip . (I5.12p admit real 
roots, we are led to study terms in the Puiseaux expansion with negative exponents. 
Consider for example the symbol 

(5.14) (1 + x n )£ m + 1, m > and n > even integers. 
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The local ellipticity condition (15. 12ft is satisfied, whereas (15. lip reduces to £ m = 0. 
We obtain 

^{x)=e j \x\- n/m + o{\x\- n/m ), J = l,...,n 

where Cj are the mth roots of — 1. Since m is even, line,- 7^ 0, and condition (ll.lOp 
is satisfied if —n/m > — 1, i.e. n < m. So we may conclude global regularity 
for the corresponding Weyl operator, in this case. When n > m, both (II. 9p and 
(ll.lOp fail. In fact, the operator is Fuchsian at 00 for m = n, and regular at 00 if 
n > m. Hence the solutions of the homogeneous equation belong to <S'(R) and do 
not belong to iS(M). We address to Camperi [3] for a general class of symbols of 
the form (I5.14p . 

Remark 5.9. According to the definition of Schwartz [17], a differential operator 
with polynomial coefficients P is hypoelliptic in M if for every open subset Q C K. 
we have: 

(5.15) u G V(Q) and Pu G C°°{Q) 4«G C°°{Q). 

We observe first that that every operator of the form (ll.ip with a mj o 7^ is obviously 
hypoelliptic. In the general case, considered in Remark I5.1| we have to take into 
account the coefficient Q(x) of the leading derivative, say of order p < m: 

Q( x ) = c p^- 

/3<m—p 

Let us write x\, . . . , x r , r < m — p, for the real zeros of Q(x). The hypoellipticity of 
P is granted in M. \ {x\, . . . , x r }. It is then easy to prove that the global regularity 
in (II. 3p implies hypoellipticity in (I5.15p . 

In fact, assume u G V'(tt) and Pu G C°°(n). Let G C °°(fi), </>(x) = 1 
in a neighborhood of the points x\,...,x r which belong to Q. We know that 
sing — suppw C {x%, . . . ,x r } D Q, hence (1 — 0)w G C°°(Q). On the other hand 
(pu G S'(Q) C S'(R) and P(<fm) = Pu - P(l - <j>)u G C^{tt) C S(R), hence 
<pu G C^°(R) by the global regularity, and therefore u G C°°(ft). 

Note that this argument cannot be extended to operators with polynomial co- 
efficients in R n , n > 1, because the manifold where the local ellipticity fails is in 
general non-compact. 

Note also that Schwartz hypoellipticity does not imply global regularity. In 
dimension 1, consider for example the operators P\ = D x , P2 = D x + x h , h > 2, 
which are obviously Schwartz hypoelliptic, but not globally regular. In particular, 
P 2 keeps the Schwartz hypoellipticity after any possible metaplectic conjugation. 
For a study of the Schwartz hypoellipticity at xi,...,x r , we address to Kannai 
|12j . where a necessary and sufficient condition was given under an asymptotic 
separation condition, similar to (II. 9p . 
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